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The  shock  Hugoniot  and  adiabatic  release  curves  of  sedimentary  rocks  display  a  range  of  different  behaviors  that 
affect  the  decay  of  shock  waves  propagating  away  from  a  confined  source.  Many  of  the  minerals  making  up  sandstone.s, 
shales,  and  limestones  undergo  phase  transitions  when  sbtxrked  to  pressures  of  interest  to  studies  of  coupling  of  energy  from 
explosive  sources  into  far  field  seismic  waves.  Both  sandstones  and  limestones  have  been  ob,served  to  exhibit  elastic 
precursors  and  multiple  wave  behavior  at  shock  vcUxiues  up  to  3,7  km/s  and  5.7  km/s,  respectively  [Ahrens  and  Cireg.son, 
1964J,  Hysteresis  in  the  shock-release  paths  of  nuienaLs  results  in  irreversible  energy  deposition,  thus  depnving  the  shock 
wave  of  energy  required  to  drive  its  propagauon  and  resulting  in  a  more  rapid  decay  and  less  efficient  coupling  to  far  field 
waves  than  expected  from  geometrical  effects  alone  Thus,  data  constraining  this  behavior  in  rocks  are  essential  for 
determining  the  cavity  volumes  for  various  degrees  of  decoupling  for  different  lithologies. 

The  behavior  of  most  sedimentary  rocks  during  adiabatic  release  has  not  been  experimentally  investigated  to  a 
significant  extent  Our  objective  was  to  obtain  expenmenial  data  for  shock  and  release  behaviors  of  elastic  and  carbonate 
sedimentary  rocks  and  use  those  data  along  with  data  from  other  sources,  to  develop  a  theoretical  model  of  the  hysiereiic 
shock-release  paths  of  these  rocks.  We  use  this  model  to  constrain  the  energy  deposited  irreversibly  in  these  rocks  by  the 
passage  of  shock  and  release  waves.  The  release  model  we  use  is  a  modification  of  that  used  by  Sekine  ei  al  [1993. 
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OBJECTIVE 


The  shock  Hugoniot  and  adiabatic  release  curves  of  seditoentarv  rocks 
display  a  range  of  different  behaviors  that  affect  the  decav  of  shock  waves 
propagating  away  form  a  confined  source.  Many  of  the  minerals  making  up 
sandstones,  shales,  and  limestones  undergo  phase  transitions  when  shocked  rc 
pressures  of  interest  to  studies  of  coupling  of  energy  from  explosive  sources 
into  far  field  seismic  waves.  Both  sandstones  and  limestones  have  been 
observed  to  exhibit  elastic  precursors  and  multiple  wave  behavior  at  shock 
velocities  up  to  3.7  km/s  and  5.7  km/s ,  respectively  [y^hrcns  and  Cregson. 
1964).  Hysteresis  in  the  shock-release  paths  of  materials  results  in 
irreversible  energy  deposition,  thus  depriving  the  shock  wave  of  energv 
required  to  drive  its  propagation  and  resulting  in  a  more  rapid  decav  and  ies 
efficient  coupling  to  far  field  waves  than  expected  from  geometrical  effects 
alone.  Thus,  data  constraining  this  behavior  in  rocks  are  essential  for 
determining  the  cavity  volumes  for  various  degrees  of  decoupling  for  differen 
lithologies . 

The  behavior  of  most  sedimentarv  rocks  during  adiabatic  release  has  not 
been  experimentally  investigated  to  a  significant  extent.  Our  objective  was 
to  obtain  experimental  data  for  shock  and  release  behaviors  of  clastic  and 
carbonate  sedimentary  rocks  and  use  those  data  along  with  data  from  other 
sources,  to  develop  a  theoretical  model  of  the  hysteretic  shock-release  paths 
of  these  rocks.  We  use  this  model  to  constrain  the  energv  deposited 
irreversibly  in  these  rocks  by  the  passage  of  shock  and  release  waves.  The 
release  model  we  use  is  a  modification  of  that  used  by  Sekine  er  al .  \19?j  . 
submiccedj  for  granite. 

KESEARCH  ACCOMPLISHED 

We  have  studied  the  shock  compression  and  adiabatic  release  behavior  of 
Coconino  sandstone  with  an  initial  density  ”  2.329  Mg/m^,  Soienhofen 
limestone  with  initial  density  Poc  “  2.613Hg/m^,  and  a  calcareous  shale  with 
initial  density  poo  -  2.574  Mg/m-.  The  sandstone  is  essentially  quartz  with 
only  minor  amounts  of  other  constituents  and  the  limestone  is  essentially  pur 
calcite  with  very  minor  amounts  of  other  constituents.  The  shale,  based  on  a 
modified  CIPW  normative  mineral  calculation,  is  composed  of  subequal  amounts 
of  quartz,  clays,  and  calcite. 

The  features  of  the  Coconino  sandstone  Hugoniot  (figure  1)  can  be 
identified  with  those  of  the  quartz  Hugoniot.  At  very  low  pressures,  there  i 
an  elastic  shock  wave  which  precedes  the  plastic  deformation  shock  wave.  The 


i 

I 


•  •••••• 


Volume  (mV^g) 


> 

Figurt  1.  Hugoniot  curve  and  adiaDatic  rexease  paths  for  Coconino  sandstone. 


limiting  stress  of  this  elastic  wave  is  the  Hugoniot  elastic  limit  (HEL). 

Above  the  HEL.  the  Hugoniot  may  be  broken  down  into  low  pressure  phase  (LPPi 
and  high  pressure  phase  (HPP)  stability  regions,  with  an  intervening  mixed  I 

phase  (MP)  region  where  insufficient  energy  is  provided  by  the  shock  wave  to 
drive  the  formation  of  the  HPP  to  completion.  Based  on  the  data  of  Ahrens  arid 
Cregson  11964],  the  HEL  of  sandstones  seem?  very  stronglv  dependent  on  the 
initial  density  of  the  rock.  The  elastic  wave  velocity,  which  has  no 
svsteaatic  dependence  on  strain  rate  at  a  given  initial  density,  varies  from 
2.8t.2km/s  for  “1.961Mg/m5  to  3.6±.lkBi/s  for  pcc  -  2  .  Hg/a- .  I 

Above  the  HEL,  the  sandstone  LPP  and  HPP  can  be  identified  with  quartz 
and  stishovite.  The  sandstone  Hugoniot  is  well  modelled  using  the  equation  of 
state  (EOSJ  parameters  of  quartz  for  the  LPP  region  and  stishovite  for  the  HPP 
region.  The  MP  region  is  modelled  by  a  volume  weighted  average  of  the  HPP  and 
LPP  Hugoniot  curves  in  the  P-V  plane; 

I 

«  fVjjip  +  (l-f)VLpp  (1) 

where  /  is  the  mass  fraction  of  HPP.  We  assume  that  f  varies  linearly  with 
pressure  from  0  at  the  lower  bound  to  the  MP  region,  where  HPP  formation 
begins,  to  1  at  the  upper  bound  of  the  MP  region,  where  HPP  formation  is 

complete.  Consideration  of  the  present  data  and  data  from  Van  Thiel  [1977]  I 

suggests  that  HPP  formation  begins  at  a  Rankine -Hugoniot  energy  of  Ppp  = 

0 . 95±.05  MJ/kg  and  is  complete  at  Frb  =  3.8±.2  MJ/kg,  independent  of  the 
initial  density  of  the  sample.  £rh  is  given  by 

£rh  -  HPbCVb-Poo) 

» 
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where  Pg  and  Vg  are  the  pressure  and  specific  volume  on  the  Hugoniot  and  V;  it 
the  initial  specific  volume  of  the  material. 

Release  adiabats  from  the  LPP  region  are  well  modelled  by  paths  which 
have  the  same  curvature  as  the  nonporous  quartz  Hugoniot  ir  the  P-T  plane.  We 
model  this  behavior  by  assuming  that  all  release  paths  from  the  LPP  region  oi 
the  Hugoniot  are  isentropes  characterized  by  the  same  zero  pressure  bulk 
modulus,  JCso  •  find  its  pressure  derivative,  Ks' .  The  zero  pressure  density  for 
fi  given  release  path  is  constrained  by  the  requirement  that  the  isentrope 
intersect  the  Hugoniot  curve  at  the  proper  point.  We  find  Murnaghan  equation 
parameters  for  these  isentropes  to  be  Xsc  '■42.7GPa  and  Xs'  -  5.0, 

Release  from  the  HPP  region  is  modelled  in  a  similar  fashion,  following 
a  curve  which  mimics  the  stishovite  region  of  the  quartz  principal  Hugoniot . 
These  isentropes  are  described  by  Murnaghan  equation  parameters  Xs:  -  577  GPa 
and  Xs'  “  3.35.  The  material  follows  these  curves  down  to  a  critical  pressure 
Pj.  where  the  HPP  begins  to  convert  to  a  diaplectic  glass  with  a  density  of 
2.27Mg/m^.  From  to  P-0,  we  assume  that  the  HPP  release  curve  follows  a 
straight  line  in  P-V  space  to  the  diaplectic  glass  volume  at  zero  pressure. 

Release  from  the  MP  region  is  modelled  as  a  mixture  of  the  LPP  and  HFF 
release  paths,  according  to  the  prescription  in  equation  (1).  with  the 
fraction  f  of  the  HPP  being  frozen  at  the  value  achieved  on  the  Hugoniot 
Based  on  the  results  of  Sekine  et  al .  [1993,  submit c ed ] .  we  take  P.  to  be 
volume-dependent : 

P;  -  32.2  -  85.51’  I  3 

where  P^  is  in  GPa  and  I’  is  the  total  specific  volume  in  units  of  my-LMg . 

The  shale  Hugoniot  curve  (figure  2)  also  has  elastic,  LPP.  ,MF .  and  HPF 
regions.  The  HEL  of  the  shale  is  -2 . 35± . 15  GPa .  corresponding  to  a  shock 
velocity  of  4.7km/s.  The  MP  region  in  the  shale,  based  on  preliminary 
results  extends  from  -15  GPa  to  -35  GPa.  Taking  a  nonporous  zero-pressure 
density  of  pc  “  2.650, Mg/m-  for  the  shale  LPP,  we  find  that  the  shock  wave  data 
can  be  fit  with  LPP  principal  isentrope  parameters  of  Xsz  -  33  6 GPa  and  - 
8.5,  under  the  assumption  that  the  thermodynamic  Gruneisen  parameter  at  (■_  is 
7-  -  1.0  and  n  -  1.  where 


y  -  yo(Po/p'>^‘  < ' 

A  preliminary  model  for  the  HPP  based  on  MP  region  data  and  the  LPP  equation 
of  state  and  assuming  the  MP  region  is  modelled  by  equation  (1).  is  p;  -  3.172 
Mg/m^,  Xso  -  70  GPa,  Xs'  -  4.0,  and  70  -1.0  with  n  -  1 .  The  energy  of 
transition,  from  low  pressure  phase  to  high  pressure  phase  at  STP  is  0.27 

MJ/kg,  based  on  the  assumption  that  the  clay  minerals  are  similar  to  muscovite 
in  that  they  do  not  undergo  a  shock-induced  phase  transition  and  that  calcite 
converts  to  the  high  pressure  phase  and  quartz  converts  to  stishovite. 

Release  paths  from  both  the  LPP  and  HPP  regions  of  the  shale  Hugoniot 
all  fall  on  or  above  the  Hugoniot,  in  contrast  with  the  sandstone  release 
paths.  We  find  that  release  paths  from  the  shale  LPP  region  are  well  modelled 
by  a  suite  of  isentropes  with  Murnaghan  equation  parameters  Xsc  -  39  GPa  and 
Xs*  -6.1  for  release  from  shock  pressures  below  14  GPa  .  Above  Pr  -  14GPa.  the 
LPP  release  paths  change  fundamentally,  being  better  described  by  a  suite  of 
isentropes  with  Xso  -  31  GPa  and  Xs'  -  1.  As  with  the  sandstone  LP'’  and  HPP.  the 
zero  pressure  density  for  each  release  path  is  constrained  by  the  requirement 
that  the  release  path  intersect  the  Hugoniot  at  the  correct  volume.  Release 
paths  from  the  HPP  region  can  be  modelled  by  a  suite  of  isentropes  with 
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Figure  2.  Hugoniot  curve  and  adiabatic  release  paths  for  calcareous  shale. 
The  large  expansions  from  the  low  pressure  phase  region  above  14  GPa  is  also 
observed  in  pure  carbonates.  The  release  paths  ab)ve  -50  GPa  are  not  well 
constrained  by  the  available  data. 


Murnaghan  parameters  ATsc  -  66  GPa  and  Ks'  -  1-  Interestingly,  the  HPP  release 
paths  do  not  show  as  rapid  expansion  as  the  LPP  release  paths  from  Pf.;  >  l4GPa. 
Similar  behavior  is  shown  by  th.  release  data  of  V’izgirda  and  Ahrens  [1982, 
for  aragonite,  suggesting  that  this  behavior  in  the  shale  can  be  attributed  to 
the  presence  of  calcite,  suggesting  that  evolution  of  CO2  may  be  important  in 
this  region.  The  MP  release  paths  are  well  modelled  by  a  mixture  of  the  HPP 
and  LPP  release  paths  according  to  equation  (1).  with  f  frozen  at  the  Hugoniot 
value . 

The  limestone  Hugoniot  is  complicated  at  low  pressures  by  an  elastic 
wave  and  several  phase  transitions,  causing  a  multiple  wave  structure  with  as 
many  as  five  separate  shock  fronts  observed  in  some  experiments  [Ahrens  and 
Cregson,  2964],  The  multiple  wave  structure  does  not  disappear  until  shock 
pressures  above  -20GPa,  corresponding  to  a  shock  wave  velocity  of  -5.9km/«. 
are  attained.  Above  20  GPa  (figure  3),  the  Hugoniot  becomes  much  less 
complicated.  A  single  high  pressure  phase  region  extends  from  31  GPa  to  over 
100  GPa.  Based  on  the  results  of  Vizgirda  and  Ahrens  {1982],  we  model  this 
phase  with  EOS  parameters  po  “  3.05  Mg/m^,  Kso  -  85.5  GPa,  -  3.75,  7c  “  143, 
n-0.6,  and  Eu  -O.lMJ/kg.  Between  10  GPa  and  31  GPa.  the  final  shock  state 
is  modelled  as  a  mixture  of  this  high  pressure  phase  and  a  lower  pressure 
phase  that  has  po  -  2 . 985  Mg/m^ ,  iCso  -  110  GPa ,  K^'  -  4,  70  -1.44,  n-0.6,  and 
Etz  “0-  These  parameters  are  not  a  true  equation  of  state,  but  provide  a 
convenient  estimate  of  the  Hugoniot  of  the  lower  pressure  phase.  As  with  the 
mixed  phase  regions  of  the  shale  and  sandstone,  we  model  this  mixture  by 
equation  (1)  with  the  fraction  f  of  the  higher  pressure  phase  varying  linearly 
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Figure  3.  Hugoniot  curve  and  adiabatic  release  paths  for  Solenhofen 
limestone.  The  shock  behavior  below  20  GPa  is  very  complicated.  The  release 
paths  are  constrained  by  the  new  data  obtained  in  this  studv  and  bv  the  data 
of  Vizgirda  and  Ahrens  12982',. 


with  pressure  between  the  phase  transition  onset  at  10  GPa  and  completion  at 
31  GPa . 

Release  paths  from  the  single  phase  region  of  the  Hugoniot  are 
tentatively  modelled  as  isentropes  with  Murnaghan  equation  parameters  Xs:  -  66 
GPa  and  Ks’  -2.3,  based  on  our  results  and  those  of  i^izgirda  and  Ahrens 
I2982j  for  aragonite.  For  the  purposes  of  calculating  the  mixed  phase  release 
paths  from  Hugoniot  states  between  14  GPa  and  31  GPa,  we  model  the  lower 
pressure  phase  release  paths  as  isentropes  with  Murnaghan  equation  parameters 
Xsc  ”  30  GPa  and  Ks'  -  1 . 

CONCLUSIONS  AND  RECOMMENDATIONS 

We  apply  the  results  for  the  hysteretic  shock-release  paths  of  the  rocks 
studied  to  calculate  the  fraction  f'  of  the  energy  ^RB  that  is  irreversibly 
deposited  in  the  material  for  a  given  shock  pressure  using  the  equation 

f  •  1  -  (£r.i/£'rh) 

where 

£r,i  -  S 

V^R 

where  the  subscripts  H  and  R  on  F  denote  the  Hugoniot  and  release  volumes, 
respectively.  The  resulting  estimates  for  f  are  presented  in  figure  4 


(5) 

(6) 


•  •••••• 


i-or 


I 


Shock  Pressure  (GPo) 


Figure  4.  Fraction  f'  of  the  Rankine-Hugoniot  energy  irreversibly  deposited 
in  sandstone,  shale,  and  limestone,  as  a  function  of  shock  pressure.  The  low- 
values  for  the  limestone  and  shale  can  apparently  be  attributed  to  the  release 
behavior  of  calcite. 


the  sandstone,  f  approaches  1  at  low  shock  pressures  because  of  pore 
collapse,  but  rapidly  drops  to  -O.hS  at  14GPa,  rising  slovly  thereafter  to 
-0.77  at  100  GPa.  Because  the  release  paths  for  both  the  shale  and  the 
limestone  show  large  expansions,  the  value  of  fgti  i®  close  to  Efi£.  As  a 
result,  f'  is  small  for  both  of  these  materials,  falling  below  0.4  for  Ph  >  3 
GPa.  At  shock  pressures  below  -30  GPa,  both  shale  and  limestone  show  regions 
where  f*  =  0,  reflecting  the  large  expansion  of  calcite  on  release  from  the 
low  pressure  region  of  the  Hugoniot .  The  shale  shows  similar  behavior  above 
70  GPa  in  the  calculations,  but  this  behavior  is  not  well  constrained  bv  the 
data  presently  in  hand.  These  results  suggest  that  shock  waves  from  sources 
confined  in  volatile-rich  rocks,  especially  those  rocks  rich  in  carbonates, 
will  decay  more  slowly  (by  a  factor  of  -3)  than  those  in  volatile-poor  rocks 
such  as  sandstones.  As  a  result  coupling  of  explosions  is  expected  to  be  jDore 
efficient  in  volatile-rich  rocks. 
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